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ABSTRACT 

We use cosmological hydrodynamical simulations including star formation and metal enrichment to 
study the evolution of the chemical properties of galaxy-like objects at high redshift in the range 
0.25 < z < 2.35 in a hierarchical clustering scenario. As the galactic objects are assembled we find 
that their gaseous components exhibit neutral Hydrogen column densities with abundances and scatter 
comparable to those observed in damped Lyman-a systems (DLAs). The unweighted mean of abun- 
dance ratios and least square linear regressions through the simulated DLAs yield intrinsic metallicity 
evolution for the [Zn/H] and [Fe/H], consistent with results obtained from similar analysis of available 
observations. Our model statistically reproduces the mild evolution detected in the metallicity of the 
neutral hydrogen content of the Universe, given by mass-weighted means, if observational constraints 
are considered (as suggested by Boissee et al. 1998). For the a-elements in the simulated DLAs, we 
find neither enhancement nor dependence on metallicity. Our results support the hypotheses that DLAs 
trace a variety of galactic objects with different formation histories and that both SNI and SNII are 
contributing to the chemical enrichment of the gas component at least since z k, 2. This study indicates 
that DLAs could be understood as the building blocks that merged to form today normal galaxies within 
a hierarchical clustering scenario. 

Subject headings: galaxies: formation- evolution - abundances - cosmology: theory - dark matter - 
methods: numerical 
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weighted by the HI column density (A^hi) ( Prochaska & 



StudioD of DLA aboorborG have provided hinto on the 



Wolfe 20001 ; |Pettini et al. 1999| ; [Vladilo et al. 2000D in con 



properties of structure at high redshift, from both kine- 
matical and chemical points of view ( |Lu et al. 199^ ; ^ei 



Itini et al. 1997|; |Wolfe et al. 1995|; Haehnelt, Steinmetz fc [Edmunds fc Phillips 1997|; |Pei et al. 1999| 



trast with some theoretical models which predict substan- 
tial evolution of the neutral hydrogen abundances (e.g., 

see Gen & Os- 



Rauch 1998; |Pettini et al. 1999| ; [Prochaska fc Wolfe 19991 ). 
Analysis of the absorbers at low z show that these sys- 



tems have diverse morphologies (e.g., Rao & Briggs 1993 



Lc Brun et al. 1997) while at high z their nature remains 



unclear. First kinematical studies suggested that high red- 
shift DLAs were large disks similar to those in present-day 
spiral galaxies (e.g., Wolfe et al. 1995; Prochaska & Wolfe 



1997). However, Haehnelt et al. (1998) show that they are 



also consistent with being the building blocks of local typi- 
cal galaxies, as predicted by hierarchical clustering scenar- 
ios (HCSs). The metallicity properties of DLAs sup port 



triker 1999 for a different result). However, the unweighted 
mean mctallicities and linear regressions through the in- 
dividual data show stronger evolution signals as reported 
in recent works (i.e., Prochaska fc Wolfe 2000; Vladilo et 
al. 2000; Hon et al. 2001). Detecting evolution in chem- 
ical abundances of HI gas mass may be difficult due to 
several biasing factors. Among them, dust depletion and 
obscuration have been probed to be a source of uncertain- 
tics for some element measurements (e.g., Vladilo 1998; 



Hou et al. 2001), while the limited redshift coverage with 



the idea that th ey are chemically young objects (e.g.. Pet 



tini et al. 1999[ , Prochaska et al. 2000). However, the de- 



pendence of their mctallicities on 2; is a controversial point 
since, until recently, only very weak or non-evolution at 
all has been reported from the analysis of the mean metal- 
licity of the neutral hydrogen with the abundance ratios 
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the known DLAs may not be adequate for detecting the 
presence of evolution. As first pointed out by Boisse et al. 
(1998), DLAs with high A'hi and high metallicity might be 
missing from the data as a consequence of dust extinction 
which would imply a biased determination of evolution. 
For the a— elements, dust-corrected data show mild en- 
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hancement or near solar values for most observed elements 
in contrast to the mctallicity pattern of metal-poor Galac- 
tic stars (e.g., [Vladilo 19981 ; [Pettini et al. 2000| ). This 
might imply a different SF history than that of the Milky 
Way, although Prochaska et al. (2000) found an agreement 
between the properties of the Galactic thick disk and the 
DLAs. 

In HCSs, where galaxies are formed by aggregation of 
smaller substructures, mergers and continuous gas infall 
play an important role, contributing to regulate the star 
formation (SF; Tissera 2000) and chemical histories (e.g., 
Cora et al. 2000) of the galactic objects. As suggested 
by Haehnelt et al. (1998) from a kinematical study of 
numerical simulations, DLAs could be the progenitor sub- 
structures of today galaxies that merge to form galactic 
ob j ec ts in a hierarchical clustering scenario. The variety 



of mc |rphological types together with the fact that they 
tend to be chemically young seem to support this hypoth- 
esis. With the aim at understanding the possible link be- 
tween the structure in a HCS and the nature of DLAs, 
in this letter, we assess the chemical properties of today- 
galaxy building blocks as a function of the redshift. For 
this purpose, we focus our study on the analysis of the 
unweighted mean metallicities of the neutral hydrogen as 
a function of z since they are more sensitive to the prop- 
erties of the individual objects. We also comment on the 
mass-weighted mean metallicities that are related to the 
chemical content of the Universe. We use hydrodynami- 
cal cosmological models that treat the non-linear evolution 
of both baryons and dark matter in a self-consistent way, 
providing a well-described history of formation, and in- 
clude star formation and chemical evolution. Results of 
the simulations are compared with available observational 
data on chemical abundances of DLAs. We discuss impli- 
cations for galaxy formation. 

2. CHEMICAL PROPERTIES OF GALACTIC OBJECTS 

Our hydrodynamical chemical simulations follow the 
joint evolution of the dark matter and baryons within a 
cosmological context (Tissera et al. 1997) including SF 
and chemical evolution. Stars are formed from cold and 
dense gas in a convergent flow according to the Schmidt 
law. Gaseous particles are gradually transformed into 
stars in different SF episodes. The contribution of type 
I (SNI) and type II (SNII) supernovas from each of these 
SF episodes to the chemical enrichment of the gas compo- 
nent are taken into account according to stellar evolution 
models and metallicity enrichment yields. Chemical ele- 
ments generated in a given particle are distributed among 
gas particles within its neighboring area, weighting each 
contribution with a kernel function that depends on the 
distance (i.e., by using the smooth particle hydrodynam- 
ics technique). 

We adopt the yields given by Woosley & Weaver (1995) 
for SNII and those given by Thielemann, Nomoto & 
Hashimoto (1993) for SNI. A time delay of 10^ yrs is as- 
sumed for binary star systems to explode as SNI. We adopt 
a fixed Salpeter Initial Mass Function with lower and up- 



per mass cut-offs of 0.1 Mq and 120 Mq, respectively. 
In this work the effects of energy injection into the inter- 
stellar medium due to supernova explosions have not been 
included. For a detailed discussion of the chemical model 
see Mosconi et al. (2000). 

We analyze cosmological simulations of a typical 10 
Mpc cube volume represented by 64'^ equal mass particles 
(rjft = 0.1). Initial conditions are consistent with a Stan- 
dard Cold Dark Matter Universe [H = 50 km s"^Mpc"^) 
with cluster abundance normalization, a% = 0.67. We have 
run a set of three simulations with different realizations 
of the power spectrum, estimating averaged results over 
them. The SN parameters adopted in these simulations 
correspond to those giving t he best agreement wi th obser- 
vations of galaxies at z = ( [Mosconi et al. 200C ) and the 



[O /F e] abundance pattern in the Milky Way (|Cora et al 



20001) 



We identified galactic objects at their virial radius at dif- 
ferent stages of evolution of the simulated volume. Galac- 
tic objects are formed by a dark matter halo and a bary- 
onic component in the form of gas and stars. To dimin- 
ish numerical resolution problems, we analyze galactic ob- 
jects with more than 200 baryonic particles within their 
virial radius and in the range 0.25 < z < 2.35. Conse- 
quently, the analyzed objects have virial velocities within 
« 100 - 250 km s^^. We study a total number of 380 
galactic objects satisfying the above conditions which are 
all considered as possible absorbers. 

It is likely that DLAs observations map the chemi- 
cal properties of the gaseous disks not expected to be 
good tracers of the metal content at the central regions 
(e.g., Jimenez et al. 1999|; ^omerville et al. 2001 



Savaglio 2000). Hence, in order to carry out a suit- 



able comparison between the simulated galactic objects 
and DLAs observations, we use a Monte Carlo technique 
to simulate line-of-sights (LOS) and estimate the chem- 
ical properties of the neutr al hydrogen compo nent with 
iVni > 2 X lO^^atoms/cm-^ ([Wolfe et al. 1986|) along the 
LOS. We assume that the hydrogen mass in a gas particle 
remains neutral if no SF activity has ever occured within 
that particle. 

As a combined result of dynamical evolution, mergers 
and interactions, the SF rate history of each galactic ob- 
ject can be described as a contribution of an ambient SF 
rate and a series of starbursts (e.g., Tissera 2000 ). The 
timings between starbursts are not ad hoc, but given nat- 
urally by the evolution of the objects in the HCS adopted. 
Thus, we have a consistent description of the chemical en- 
richment of the stellar populations and gaseous component 
due to the fact that the different ejecta times of SNI and 
SNII can be properly taken into account. Note that the 
simulated sample is not affected by dust depletion or ob- 
scuration. These facts make our simulations powerful tools 
to explore the nature of DLAs. 

For each galactic object selected at a given redshift we 
define the HI mass-weighted mean of abundance ratio of 
elements K and J, [K/J] in the column of HI along a LOS, 
as: 
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[K/J] = log 



log(K/J)e 



(1) 



where rip is the total number of gas particles belonging 
to a galactic object along a certain LOS, Mf^ the neu- 
tral hydrogen mass of the z*'* particle, K, and their 
chemical abundances, and (K/J)q the corresponding solar 
abundance ratio. Similarly, we can define the mean abun- 
dance ratios for the stellar populations. This point will be 
discussed in a separate paper. 

In order to study the metallicity evolution of DLAs, we 
consider the abundance ratios of [Zn/H] and [Fe/H]. In our 
sample both elements can be measured without introduc- 
ing biases such as those produced by dust extinction and 
depletion, or, in the case of zinc by the fact that this ele- 
ment is difficult to detect at low metallicities. Note that 
the zinc is usually considered a reliable tracer of metallicity 
because it is weakly depleted onto dust grains, h owever, its 
nucle osynthesi s remains to be fully unders tood ( Prochaska 
& Wc |lfc 2000| ; iPrantzos fc Bossier 200"o|) . Current mod- 



els for the production of this element ( Woosley fc Weaver 
19950 have problem in reproducing the observational fact 
that, regardless of metallicity, stars in the halo and thin 
disk of the Galaxy have nearly solar [Zn/Fe] values (see 
also Hou et al. 2001). Nevertheless, recent results sug- 
gest that Zn is enha nced relative to Fe in th e stars of the 
Galactic thick disk ( Prochaska et al. 2000|) These make 
results from Zn analysis more difficult to interpret. For 
the purpose of carrying out a reliable comparison between 
observations and simulations, we gathered the available 
observations from the literature (see Figure 1) and esti- 
mated the same parameters calculated for the simulations 
within the same redshift range. 



In Figure la, we show the unweighted mean 
< [Zn/HJ > = l/n-^„ [Zn/HJ ratios (where n is the num- 
ber of galactic objects in the total simulated sample that 
are located at a certain z) for the HI component along 
the LOS of galactic objects in the three simulations as 
a function of redshift. Error bars correspond to the dis- 
persion of the ratios of the galactic objects at each an- 
alyzed z. As it can be seen, the simulated neutral gas 
abundance ratios are within the observed range for DLAs. 
In order to quantify their evolution with z, we compute 
the unweighted mean < [Zn/H] > in two redshift inter- 
vals, ziow = [0.26,1.5] and Zinter = (1.5,2.35]. We ob- 
tain < [Zn/H]iow > = -0.53 ±0.06 and < [Zn/H]intor > = 
— 1.12 ± 0.17 ^. These results imply an evolution of 
0.58 ± 0.13 dex between the low and intermediate red- 
shift bins in the models, while DLAs observations show a 
variation of 0.22 ± 0.09 for the same redshift bins. The 
observed mean values are < [Zn/H]iow > = —0.89 ± 0.10 
and < [Zn/H]intCT > = -l.U ± 0.07. Hence, at high z the 
simulated and observational values agree very well. How- 
ever, the low-z values differ in « 0.40 dex. This difference 
in the lower z interval could be attributed to the obser- 
vational bias first discussed by Boisse et al. (1998) that 

^Standard deviations have been calculated by using the resampling 



would prevent high iVni systems with high metallicities 
to be detected (see also Prantzos fc Bossier 2000| ; Savaglio 

2000) . By analyzing the chemical properties of the individ- 
ual DLAs, Vladilo et al. (2000) found an anticorrelation 
signal for the [Zn/H] data which yields a linear regression 
with a slope of dlogZn/dz ~ —0.32 ± 0.13, detecting in- 
trinsic evolution for z < 3.5. A similar analysis applied 
to our set of observational data within our z-range yields 
dlogZn/dz = -0.24 ± 0.11. For the simulated DLAs, we 
obtained a slope of dlogZn/dz = — 0.42±0.08 which, albeit 
larger, is statistically consistent with observations. 

Figure lb shows the unweighted < [Fe/H] > versus z for 
the same galactic objects plotted in Figure la. This ele- 
ment suffers strong gas depletion but can be observed over 
a larger z interval, its nucleosynthesis is better understood 
and can be easily detected at low metallicities. ^From Fig- 
urelb we see that the estimated ratios for the HI along the 
LOS in the simulations match the observational range of 
DLAs. Estimates of their unweighted < [Fe/H] > at ziow 
and Zintcr give a 0.67 ± 0.13 dex evolution in the simula- 
tions. Calculations for the available DLAs observations in 
the range 0.25 < z < 2.35 show, in this case, a signal for 
evolution of 0.48 ± 0.13 dex. A linear regression through 
the observed [Fe/H] data within the same z-range shows 
a slope of —0.50 ± 0.17 (see also Savaglio 2000; Hou et al. 

2001) . The linear correlation of the simulated [Fe/H] ra- 
tios gives dlogFe/dz — —0.50 ± 0.09. All these values are 
in suitable agreement within the statistical uncertainties. 

The mass-weighted mean metallicities of HI assess the 
evolution of the chemical properties of the population as 
a whole and, consequently, of the chemical history of the 
Universe. As previously reported, both the observed A^hi- 
weighted mean [Zn/H] and [Fe/H] ratios show mild or no 
evolution for z < 3.5 (Pettini et al. 1999; Prochaska & 



Wolfe 2000), while most theoretical models predict sub- 



stantial change of the metallicity of neutral gas with red- 
shift. Estimates of the corresponding ratios in our model 
at low and intermediate z intervals also yield a signal of 
evolution: 0.33 ± 0.08 and 0.47 ± 0.09 for the [Zn/H] and 
[Fe/H], respectively. 

We would like to stress that the conclusions about the 
evolution of the observed [Zn/H] and [Fe/H] by using mean 
values at low and intermediate z should be taken with 
caution because of the small number statistics in the low 
z-bin. More observations of DLAs at low z are needed 
to draw final conclusion on this point. But, mainly, more 
observations of DLAs with high iVni and high metallicity 
are required to properly assess the presence of evolution. 
In fact, if DLAs with these characteristics were missing 
from the data then the presence of evolution could be dif- 
ficult to detect. In order to mimic this effect, Prantzos & 
Boissier (2000) applied a filter to their model by select- 
ing the DLAs according to the observationally determined 
constraint 18.8 < [Zn/H] -f logA^ni < 21. Following their 
work, we estimate the evolution of the filtered simulated 
DLAs, using the [Zn/H] ratios. For the unweighted means 

bootstrap technique with 500 random samples. 
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we obtained an evolution signal of 0.46 ± 0.10 while the 
weighted ratios show a weaker change consistent with mild 
evolution: 0.24 ±0.12. Clearly, selection biases could have 
strong implications when estimating changes with redshift. 

The standard dispersions shown in Figure 1 for both el- 
ements in the simulated sample reproduce fairly well the 
observed dispersions of DLAs: for the [Zn/H] the building 
blocks show a dispersion of (Tz,^^ — 0.68 and crzi„tor — 0-78 
while observations have Cziow — 0.42 and crzi„tor — 0.86. 
The simulated [Fe/H] ratios have Czj^^ = 0.60 and crzi„t„, — 



0.77 compared to the observed ones of CTzj 



0.45 and 



"'zintor = 1.24. In our models this dispersion arises as the 
results of the variety of galactic objects with diverse evo- 
lutionary histories: merger tree, star formation, chemical 
enrichment, etc., at different z, supporting a similar origin 
for the large dispersion observed in DLAs. 

We estimate the unweighted mean [Si/Fe] and [S/Fe] ra- 
tios in order to analyze the behavior of the a-elements in 
the HI gas mass of the simulated galactic objects. ^From 
Figure 2a we can appreciate that the mean [Si/Fe] ratios 
have nearly solar values and show no dependence on metal- 
licity, consistent with the dust-corrected observations re- 
ported by Vladilo (1998) and Pettini et al. (2000) who 
assumed that Zn is undepleted and traces Fe. A similar 
behavior was found for mean [S/Fe] ratios. 

In Figure 2b we show the unweighted < [Zn/Fe] > abun- 
dances versus < [Fe/H] > . We find that the HI compo- 
nent of the simulated DLAs tend to exhibit [Zn/Fe] values 
around a mean of 0.27 ± 0.28 (Icr deviation). It should be 
noted that these values are consistent with the observed 
range for DLAs with no dust corrections (with a mean 
of 0.47 ± 0.31) while dust-corrected ratios show nearly 
solar values. As more extensively discussed by Prantzos 



& Boissier (2000; see also Lu et al. 1996 and Prochaska 
& Wo lfe 20"00| ) current nucleosynthesis models of the zinc 
have problem in reproducing the observed abundance pat- 
tern of the Milky- Way. If we had assumed th at Zn traces 
Fe for all metallicites (e.g., Pettini et al. 2000 ), our models 
would have produced mean ratios in agreement with those 
of Vladilo (1998) corrected by dust-depletion. 

3. DISCUSSION AND CONCLUSIONS 

Observational studies of the chemical properties of 
DLAs have provided unvaluable information on galaxy 
evolution up to z « 4. This letter aims to provide a 
theoretical frame to understand these observations within 
HCSs. The chemical hydrodynamical cosmological models 
analyzed take into account different physical mechanisms 
controlling the formation and evolution of galaxies: grav- 
itatio nal collapse, tidal torques, mergers, interactions, gas 



inf all , I gradient pressure, radiative cooling. SF and chem- 



ical enrichment by SNI and SNII explosions. In this sce- 
nario we follow the formation and evolution of galaxy-like 
objects, analyzing the chemical properties of the different 
substructures that merge hierarchically. 

In order to properly compare the galactic objects in the 
simulations with the observed DLAs, we estimated the 
abundances of HI along LOS randomly distributed. We 



found that the simulated galactic objects have gaseous 
components that, when randomly sampled, have [Zn/H] 
and [Fe/H] abundances within the observed range for 
DLAs. In the simulations, the unweighted [Zn/H] and 
[Fe/H] ratios show evolution with z in statistical agree- 
ment with results obtained from similar analysis of the 
available observations (see also Vladilo et al. 2000; 
Prochaska fc Wolfe 200"o|; Hou et al. 2001). However, 



the mass-weighted metallicities of the simulated and ob- 
served Nm show opposite trends; while the former show 
evolution, the later do not. Part of this discrepancy can be 
attributed to observational biases that might prevent low- 
density metal-poor and very metal-rich and dense zones to 
be detected by current observations. In fact, taking into 
account these effects by applying a filter to the simulated 
[Zn/H], we found, on one hand, that the unweighted means 
still show a clear evolution, and, on the other, that the 
filtered weighted abundances yield a weaker dependence 
with z in agreement with corresponding observations (e.g., 
Pettini et al. 1999). Accordingly to our results, this obser- 
vational bias seems to affect more strongly the weighted 
means suggesting that the interpretation of the evolution 
of the chemical content of the Universe by using DLAs 
should be taken with extra caution. More observations, 
principally of high density and high metallicity HI column 
densities, are needed to draw a final conclusion. 

The lack of a— enhancements in the mean abundances of 
simulated galactic objects shows that, in the redshift range 
analyzed, there are already significant contributions from 
both types of SNs, so that the mean unweighted metallic- 
ities of the iVni are nearly solar, in agreement with DLAs 
observations after dust corrections. We detect no evolution 
of [a/Fe] with metallicity implying that, for our galactic 
objects, their SF histories are the result of the superpo- 
sition of different starbursts occuring at different epochs. 
Moreover, the starbursts could have been triggered in dif- 
ferent substructures that merge to form the galactic ob- 
jects analyzed at a given z. In the case of Zn we found 
enhancement with respect to the Fe but no dependence on 
metallicity. A lower Zn nucleosynthesis production in the 
models (or assuming that Zn traces Fe) would have pro- 
duced [Zn/Fe] va lues in agreem ent with the dust-corrected 
observed DLAs ( Vladilo 1998 ). This lower Zn production 
does not affect the agreement found between the observed 
and simulated [Zn/H] at low z shown in Figure la. At high 
z, models and observations can be reconciled if one takes 
into account that several high z estimations of [Zn/H] in 
observed DLAs are upper limits instead of accurate values. 
The zinc remains a controversial element as pointed out by 
several authors ( Lu et al. 1996|: 



tini et al. 200C 



Pettini et al. 1999; Pet- 



Prantzos fc Bossier 2000; Savaglio 2000) 



whose production and evolution are not fully understood. 

To sum up, the chemical properties of galactic objects in 
our models at a given z are the result of their past evolu- 
tion: collapse, mergers, interactions, SF history, etc. Since 
we do not impose any particular constraint to select the 
galactic objects in the models, the results of our chemical 
analysis seem to support that DLAs could be the build- 
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ing blocks of today normal galaxies within a hierarchical 
context. 
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Fig. 1. — [Zn/H] (a) and [Fe/H] (b) unweighted mean abundances for the neutral hydrogen along LOS in the galactic objects in the three 
simulations as a function of the redshift (filled circles). Error bars correspond to la standard deviation. We include observational data (open 
pentagons) for DLAs taken from Lu et al. (1996), Pettini et al. (1997; 1999; 2000), Vladilo (1998) and Prochaska & Wolfe (1999; 2000). The 
downward- (upward-) pointing arrows indicate upper (lower) limits of [Zn/H]. Unweighted mean values for the simulated DLAs (dotted-open 
circles) and the observations (dotted-open pentagons) in the low and intermediate redshift intervals are also plotted. These symbols are 
superimposed in the intermediate redshift interval in Figure (a). Solid lines are the least square linear regression for the whole sample of 
simulated DLAs. 



Fig. 2. — [Si/Fe] (a) and [Zn/Fe] (b) unweighted mean abundances for neutral hydrogen along LOS in the simulated DLAs in the three 
simulations for each redshift analyzed as a function of the metallicity [Fe/H] (filled circles). Error bars correspond to la standcird deviation. 
We include observational data for DLAs without dust corrections (open pentagons) taken from Lu et al. (1996) and Prochaska & Wolfe 
(1999); dust corrected values (dotted-open pentagons) are obtained from Vladilo (1998) and Pettini et al. (2000). The downward- (upward-) 
pointing arrows indicate upper (lower) limits of [Si/Fe] and [Zn/Fe]. 



